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ABSTRACT 

We use high-resolution cosmological zoom simulations with ~ 200 pc resolution at 
z — 2 and various prescriptions for galactic outflows in order to explore the impact of 
winds on the morphological, dynamical, and structural properties of eight individual 
z = 2 galaxies. We present a detailed comparison to spatially- and spectrally-resolved 
Ha and other observations of z ~ 2 galaxies. We find that simulations without winds 
produce massive, compact galaxies with low gas fractions, super-solar metallicities, 
high bulge fractions, and much of the star formation concentrated within the inner 
kpc. Strong winds are required to maintain high gas fractions, redistribute star-forming 
gas over larger scales, and increase the velocity dispersion of simulated galaxies, more 
in agreement with the large, extended, turbulent disks typical of high-redshift star- 
forming galaxies. Winds also suppress early star formation to produce high-redshift 
cosmic star formation efficiencies in better agreement with observations. Sizes, ro- 
tation velocities, and velocity dispersions all scale with stellar mass in accord with 
observations. We find a diversity of morphological characteristics - among our two 
most massive galaxies, one resembles a quiescent grand-design spiral while the other 
is a clumpy disk undergoing a minor merger; the clumps are evident in Ha but not in 
the stars. Some galaxies have stellar densities above the threshold for compact ellipti- 
cals seen at these redshifts, although these galaxies are star- forming since we do not 
include quenching feedback. Rotation curves are generally slowly rising, particularly 
when calculated using azimuthal velocities rather than enclosed mass. Our results are 
broadly though not optimally resolution-converged. These results show that cosmolog- 
ical simulations including outflows can produce disk galaxies similar to those observed 
during the peak epoch of cosmic galaxy growth. 



1 INTRODUCTION 



The epoch around redshift z ~ 2 is the most ac- 
tive period of cosmic st ar formation (|Madau et al.l |l996; 
iHopkins &; Beacoml l2006h , and appears to be the period 
when the familiar Hubble sequence first began to emerge. 
It is thus a critical epoch for understanding how galaxies 
form, grow, and evolve into the populations we see today. 
The advent of near-infrared integral field spectrometers on 
8-10 m class ground-based telescopes has enabled spatially 
and spectrally resolved observati ons of an increasing number 
of galaxies at this epoch (e.g., [F orster Schrei ber et al.ll 2006 
Cresci et all I2QQ9I: Forster Schreiber et al. 2009; Law et al. 



2009 : IWright et al.l 120091 : lAlaghband-Zadeh et al.l 12012, 
Swinbank et al.l l2Q12h . These studies are complemented by 



space-based observations with Hubble, enabling a detailed 



study of the distributions of star formation and stellar pop- 
ulations, and Spitzer and Herschel, providing infrared data 
that better constrains the total stellar mass and bolometric 
emission of these galaxies. Such data are providing a compre- 
hensive view of the structure, kinematics, and star formation 
properties of galaxies at early stages of their evolution. 

These observations provide new and unexplored avenues 
with which to constrain models of galaxy formation and evo- 
lution. In particular, the resolved (both spatially and spec- 
trally) information that is now available provides detailed 
constraints on the assembly of galaxies at z ~ 2. The galaxy 
population at this epoch displays a remarkable level of di- 
versity, with many properties unlike anything seen locally, 
from extremely compact and (relatively ) quies cent ellipticals 
(|van Dokkum et al.1 [20081 : iBarro et ai]l2012h . to extended 



© 0000 RAS 



2 D. Angles- Alcazar, R. Dave, F. Ozel & B. Oppenheimer 



turbulent and clu mpy disks (|Forster Schreiber et al.l 12009; 
IGenzel et al.ll201lh , to galaxies formi ng stars at hundreds to 
thousands of solar masses per year (IChapman et al.ll2010l : 
lAlaghband-Zadeh et al"]l2012l : ITargett et al.ll2Q12h . Captur- 
ing this diversity, both qualitatively and quantitatively, 
presents an enormous challenge for galaxy formation mod- 
els, one that is only now beginning to be addressed using 
the latest computational tools. 

Cosmological hydrodynamic simulations have now ma- 
tured to the point that they yield galaxy populations 
broadly in agree ment with observat ions across a range 
of redshifts (e.g., IPave et~aL I l2011al lbh. While there are 
still many discrepanci es with even the latest models (e.g., 
IWeinmann et al.l 12012). overall it appears that the galaxy 
population is in general agreement with expectations from 
a A cold dark matter (ACDM) cosmology, combined with 
numerous feedback mechanisms from star formation, photo- 
ionising radiation, and active galactic nuclei (AGN) that 
help establish the properties of galaxies assembling within 
hierarchically- growing halos. The emerging paradigm is that 
galaxy evolution at this epoch is governed by a balance 
between inflows from the intergalactic medium (IGM) and 
powerful, ubiquitous outflows that intimately connect galax- 
ies and their surroundings in a " baryon cycle" of exchanging 
mass, energy, and metals fe.g. JPave et al.ll2012h . 

Resolved studies of distant galaxies present a new chal- 
lenge to such models. In order to meet this challenge, nu- 
mericists have begun to employ the "zoom" technique to 
expand the dynamic range sufficiently to model such obser- 
vations. In zoom simulations, a sub- volume extracted from 
a larger volume is re-simulated at significantly higher res- 
olution, providing a substantial increase in resolution at a 
manageable computational cost. Such simulations can simul- 
taneously capture the crucial baryon cycle on larger scales, 
while still achieving sufficiently high enough resolution to ro- 
bustly model the internal structure and dynamics of galax- 
ies. 

While inflows are generally well-predicted within the 
ACDM paradigm (albeit difficult to detect), galactic out- 
flows remain a highly uncertain and poorly-constrained in- 
gredient in galaxy formation models. It is now evident 
that powerful galactic outflows are ubiquitous in high red- 
shift galaxies, as indicated by the high frequ ency of blue- 
shifted rest-frame UV absorption lines (e.g. , 
2009; ISteidel et al] l2010l : iKornei et al] [2012; 



Weiner et al 



2011 



2012) and broad Ha emissio n-line profiles (e.g. 



Martin et al 



Genzel et al 



_ / I Newman et alj l2012h . These wind s carry out masses 
comparable to the star formation rates ([Steidel et al.ll201ol ; 
IGenzel et alboilh . though these estimates can be uncertain 
by an order of magnitude. Hence outflows are likely to play 
a central role in the evolution of galaxies and the IGM. 

Theoretically, outflows are the primary candidate for 
regulating the baryon and metal content of galaxies, while 
concurrently explaining the enrichment of the IGM. The 
physical origin and launching mechanism of such out- 
flows is still in debate, and it is usually attributed to 
energy and/or momentum input from supe rnovae (SNe) 
and/or rad iation pressure from massive stars (Mu rray et al.l 
2005. 2010|V Galaxy scale an d cosmological simulations by 
Springel fe HernauistI (|2QQ3al lbh showed that the injection of 
energy from SNe in the form of kinetic outflows with con- 
stant velocity can be effective in removing gas from galaxies, 



potentially solving the overcooling problem by regulating 
star formation rates (SFRs) and enriching the IGM. 

However, observations of dwarf starbursts and low red- 



shift luminous infrared galaxies ([Martini 20051; IRupke et all 
2005) as well as higher redshift galaxies (jweiner et al.ll20Q9h 
suggest that the properties of galactic outflows scale with 
galaxy properties - galaxies with higher masses and SFRs 
drive faster and more energetic winds - in broad agree ment 
with predictions of momentum- driven models (Mu rray et al.l 
2005). In the momentum-driven wind scenario, radiation 
from massive stars is absorbed by dust that collisionally cou- 
ples to the gas, resulting in galactic outflows for which the 
velocity and mass loading factor (i.e., the mass loss rate rel- 
ative to the SFR) scale linearly and inversely, respectively, 
with the circular velocity of galaxies. 

iQppenheimer fe Pavel {2006, 2008) implemented a vari- 
ety of outflow models into cosmological hydrodynamic sim- 
ulations of galaxy formation, and found that the scalings 
arising for momentum-driven winds yielded a signifi c ant im - 
provement over the original Springel k Her nauistI (2003b) 
"constant wind" models towards matching a wide range of 
observables, including (i) the chemical enrichment of the 



IGM at z > 2 (Oppenheimer k Pavel 2006. 


2008h. (ii) the lu- 


minosity function of high-redshift galaxies ( 


Pave et al.ll2006l; 



the galaxy 



mass- m et allicity relation ([Finlator &; Pav e 2008; P ave et al.l 
l2011bh . However, with spatial resolution of typically several 
kpc, these simulations were unable to examine the internal 
structural properties of galaxies on sub-kpc scales. 



In this work, we use cosmological hydrodynamic zoom 
simulations includi ng galactic outflows using the exact same 
prescriptions as in IPave et al.l (|2011ah to make predictions 
for i nternal structure and dyna mics of 8 re-simulated galax- 
ies ([Angles- Alcazar et a l. 2013). Our modeling is similar to 
that used in iGenel et al.l (|2012h , who f ocused on studying 
the properties of star-formin g clumps ([Genzel et al ] l201ll : 
IForster Schreiber et al.ll201lh in models with and without 
winds, finding that a similar wind prescription provided 
both realistic suppression of star formation while disrupt- 
ing clumps on short timescales. Here, we examine the im- 
pact of galactic outflows on the morphological, dynamical, 
and star formation properties in a sample of 8 re-simulated 
central galaxies in two re- simulated regions, and compare 
these properties to observations from various surveys, fo- 
cusing particularly on the Spectroscopic Imaging survey in 
the Near-infrared with SINFO NI (SINS) of z ~ 2 galaxies 
([Forster Schreiber et aLll2Q09h . 



We begin by describing our simulations in section [2] and 
present and overview of the sample of simulated galaxies in 
section [3] We analyze the impact of galactic outflows on 
the time evolution and radial structure of three example 
galaxies in section 2] where we also evaluate the effects of 
different wind models on global properties of galaxies such 
as halo baryonic fractions and cosmic star formation effi- 
ciencies, and their rotation curves. In section [5] we compare 
the properties of our simulated galaxies to available spa- 
tially and spectrally resolved observations of z ~ 2 galaxies. 
We present resolution convergence tests of our key results in 
section [6] and we summarize our results in section [7] 
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2 SIMULATIONS 
2.1 Simulation code 

Our simulations were run with an extended version of the 
N-body + smoothed particle hydrodynamic s (SPH) cosmo - 
logical galaxy formation code Gadget- 2 (!SpringeJl2005h . 
Gad get- 2 combines an entropy -conserving formulation of 
SPH (Spr ingel fe Herna uist 2002) along with a tree-particle- 
mesh algorithm for computing gravitational forces. Addi- 
tions to the public version of t he code have been described 
in lO ppenheimer fe Pavel |2008) and include models for gas 
cooling, star- format ion, chemical enrichment, and galactic 
outflows. 

We include photoionization heating starting at z = 
9 via a spatially unifor m, optically thin UV background 
(|Haardt fe Madaull200lh . We account for radiative cooling 
from primo r dial g as assuming ionization equilibrium as in 
iKatz et alJ (| 19961 ) and metal-line co oling using the colli- 
sional ionization equilibrium tables of ISutherland fe Dopital 
(|l993l ). We track the production of four metal species (C, 
O, Si, and Fe) from Type II SNe, Type la SNe, and asymp- 
totic giant branch (AGB) s tars using met all i city- d ependent 
yields as described in Qppe nheimer fe Pavel (|2008l ) . We also 
account for energy feedback from Type II and Type la SNe, 
and mass-loss from AGB stars. 

St ar formation is model e d usin g the sub-grid prescrip- 
tion of lSpringel fe Hernquistl ([20 03a). Gas particles that are 
sufficiently dense to become Jeans unstable (n > 0.13 cm -3 ) 
are treated as a two-phase interstellar medium (ISM) con- 
sisting of hot gas that condens es into cold star-formin g 
clouds via a thermal instability (Mc Kee fe Ostriked[l977h . 
Stars form from the cold phase and thermal feedback from 
Type II SNe causes the evaporation of cold clumps into the 
hot medium. Star formation is implemented probabilisti- 
cally, such that at any given step, a sufficiently dens e gas 
particle can spawn a star based on a ISchmidtl (|l959h law. 
The star particle's mass is half the original gas particle mass. 
The re sulting SFRs are t uned to be in accord with the ob- 
served [Kennicutt] (|l998aT ) relation. 

Our simulations include a galactic outflow mechanism 
that imparts kinetic energy to gas particles. The outflow is 
described by two parameters, the wind speed v w and the 
mass loading factor 77, which is the mass outflow rate in 
units of the SFR. If a particle is eligible to form stars, it is 
likewise eligible to be kicked into an outflow, with a prob- 
ability given by 77 times the star formation probability. If 
selected, the gas particle is kicked with a velocity v w in the 
direction vxa, where v and a are the particle's instanta- 
neous velocity and acceleration. Hydrodynamic forces are 
turned off until the particle reaches a density of 0.013 cm -3 , 
or a time of 20 kpc/v w has passed. This outflow implemen- 
tation i s identical to that used in larger-scale simulations 
Dave et al Tj2011al), and is very similar to that in 



of. 



ISoringel fe Hernquistl (|2003bh . 

Choices for v w and 77 thus define our "wind model". 
Throughout this paper we compare the following wind mod- 
els: 

• No wind model (nw): We turn-off galactic winds (77 = 
0). 

• Constant wind model (cw): v w — 680 km s 1 and 77 = 2 
for all galaxies. 



• Momentum- driven wind model (vzw): The kick veloc- 
ity scales with galactic velocity dispersion, a, and the mass 
loadi ng factor scales as 1/ cr, as in the momentum-conserving 
case jMurrav et al.1 2005). Following Qppenh eimer fe Pavel 
( 2008) we take v w = 3a V/l — 1 and 77 = ao/cr, where /l 
is the (metallicity dependent) critical luminosity necessary 
to expel gas from the galaxy and o~o = 150kms _1 . We run 
an on-the-fly galaxy finder to calculate galaxy m asses that 
are t hen converted to a using standard relations (Mo et all 



2.2 Simulation Runs and Analysis 

Current high resolution near-infrared observations being 
carried out by adaptive optics-assisted 10 m-class telescopes 
as well as the Hubble Space Telescope are able to resolve 
~ 1 kpc scales at z ~ 2 (e.g.. iForster Schreiber et~aDl201ll : 
iKartaltepe et al.l l2012h . In order to get full advantage of 
these observations, we used the "zoom-in" technique (e.g., 
iNavarro fe White! 1 19941 ) to carry out cosmological simula- 
tions that follow the evolution of galaxies down to z ~ 2 
with sub-kiloparsec spatial r esolution. The simulations p re- 
sented here have been used in lAngles- Alcazar et al.l (|2013l ) in 
a different context, where we constrain the growth of massive 
black holes at the centers of galaxies on cosmological time 
scales and discuss the implications of the observed black 
hole-galaxy correlations. 

We selected sizable regions for re-simulation surround- 
ing two galaxies characterized by similar stellar masses 
but different morphologies and merger histories from an 
intermediate-resolution full cosmological simulation (includ- 
ing momentum-driven winds) with 2 x 256 3 gas+dark 
matter particles in a [24/i -1 Mpc1 3 box. As described in 
Angles- Alcazar et a zoom initial conditions were 

generated by identifying all dark matter particles within the 
virial radius of each selected z = 2 galaxy, tracing their posi- 
tions back to their locations on the initial grid. Refinement 
regions were initially defined by the grid cells containing 
these particles and, subsequently, significantly enlarged by 
an iterative cleaning procedure incorporating an increasing 
number of neighboring cells. High resolution regions were 
populated with a large number of lower mass particles (yield- 
ing x64 mass resolution increase for our highest resolution 
simulations) and small scale power spectrum fluctuations 
were applied according to the spatial resolution of the re- 
fined grid. Additionally, two nested concentric layers of pro- 
gressively lower resolution were defined surrounding the high 
resolution regions in order to reduce numerical artifacts due 
to the difference in particle masses and to ensure that the 
large scale gravitational torques acting on target halos are 
accurately represented. 

The sample of galaxies presented here consists of two 
galaxies located at the centers of each re-simulated region to- 
gether with 6 additional galaxies that were found well within 
the high resolution regions and with no contamination of low 
resolution particles within their virial radius. Even though 
our sample is not mass nor volume complete, our galax- 
ies are representative of "normal" z — 2 systems. Our 
zoom simulations have (high resolution) gas particle mass 
mgas ~ 2.3 x 1O 5 M0, dark matter particle mass todm ~ 
1.2 x 10 6 M©, and softening length e « 0.47 h~ x kpc comov- 
ing (~ 224 pc physical at z = 2), equivalent to 2 x 1024 3 par- 
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Figure 1. Star formation rate surface density maps color-coded according to the projected density- weighted temperature for our 8 re- 
simulated galaxies at z = 2. For each galaxy, we show the projected gas distributions corresponding to simulations including momentum- 
driven winds (top), constant winds (middle), and no winds (bottom). Face-on and edge-on views of each galaxy correspond to the total 
angular momentum of the gas component. Galaxies are ordered from left to right for increasing stellar mass according to simulations 
with momentum-driven winds (with total halo masses ranging from ~ 10 11 to 2 x 10 12 M ). Total SFRs are listed for all galaxies and 
wind models. The region shown is 30 kpc (physical) across. 



tides homogeneously distributed in a [24/i _1 Mpc] 3 box. Ad- 
ditionally, we run simulations with a factor of 2 lower spatial 
resolution and a factor of 8 lower mass resolution in order to 
test our results for n umerical converg ence. Throughout this 
paper we assume a IChabrierl {2003) initial mass function 
(IMF) and a ACDM concordance cosmology with parame- 
ters Ov 0.72, Q M 0.28, Q h 0.046, h 0.7, a 8 0.82, 
and n = 0.96, consistent with five-year WMAP data com- 
bined with baryon aco ustic oscillations and SNe constraints 
(Komatsu et al. 2009) as well as the final nine- year WMAP 
data (|Hinshaw et al.ll2012h . 

We wrote snapshot files at time intervals ranging from 
~ 5 to 25 Myr during the simulation. Galaxies were iden- 
tified by means of the Spline Kernel Interpolative Denmax 
algorithm (skieQ) at all available snapshots independently, 
and are thus defined as bound groups of star-forming gas 
particles (i.e., gas particles with densities abov e the thresh- 
old f or star formation) and star particles (see iKeres et al.l 
2005). We used a spherical overdensity algorithm to asso- 
ciate each SKiD-identified galaxy with a dark matter halo, 
where the viri al radius was denned to enclose a mean den- 
sity given by iKitavama &; Sutol (| 19961 ) . Overlapping halos 
were grouped together so that every final halo has a central 
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galaxy and a number of satellite galaxies by construction. 
The 8 galaxies analyzed here are centrals. 

The full evolution of each z = 2 galaxy was recon- 
structed by identifying its most massive progenitor at all 
available redshift snapshots. Galaxies were also identified 
across simulations with different wind models to allow for a 
detail ed model comparison. Following I Angles- Alcazar et al.l 
(|2013r ) , we calculated structural and kinematic properties of 
each galaxy relative to the position of its most bound gas 
particle, which is a more meaningful and more stable defi- 
nition of the nominal center of the galaxy compared to that 
computed by SKID, especially during close galaxy encounters 
and galaxy mergers. 



3 GALAXY SAMPLE 

We begin by describing the sample of 8 galaxies used in 
this work. At z = 2, the masses of dark matter halos range 
from ~ 10 11 to 2 x 10 12 M© and are typically resolved with 
~ 2 x 10 5 to 4 x 10 6 particles. As we describe in detail 
throughout this paper, galaxy properties are highly depen- 
dent on the model adopted for galactic outflows. For our 
fiducial simulations including momentum-driven winds (vzw 
model), galaxies have stellar masses ranging from ~ 1.6 x 10 9 
to 3.3 x 10 10 M and SFRs in the range - 1-40 M yr" 1 at 
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Figure 2. Same as Figure [T] for the projected stellar surface density. Total stellar masses are listed for all galaxies and wind models. 



z — 2. Figure \T\ shows face-on and edge-on views of the 
star- forming gas of galaxies at z = 2. Our galaxy sample 
spans a wide range of sizes and morphologies, with the star- 
forming gas extending up to scales of 2-15 kpc from their 
centers. Some galaxies are located in high density environ- 
ments and undergo frequent interactions and mergers (e.g., 
galaxy g222) while others are characterized by a smoother 
evolution, in lower density environments (e.g., galaxy g54). 
Despite this, most galaxies appear to be rotationally sup- 
ported disks at z = 2. Figure [2] shows the projected stel- 
lar distribution corresponding to the same galaxies. Stellar 
disks can also be visually identified for most galaxies but in 
this case a higher fraction of stars seems to correspond to a 
spheroidal component. 

Galaxies are identified across the various wind simula- 
tions via their halos, making it possible to analyze the ef- 
fects of galactic winds in a galaxy- by- galaxy basis. However, 
when comparing the morphologies of galaxies across wind 
models it is important to note that galaxy interactions and 
mergers are inherently random processes. Small deviations 
of orbital parameters in different simulations can result in 
rather different morphologies at a given time, and therefore 
it is not trivial to compare morphological features such as 
spiral arms or tidal tails among the different models (e.g., 
galaxies g2743, g85, and g54). Despite this, a quick visual 
inspection of Figures [1] and [2] reveals significant differences 
between wind prescriptions. Constant winds are very effi- 
cient in removing gas from galaxies and result in lower mass 
galaxies with typically lower SFR surface densities compared 
to the momentum-driven winds. This is particularly evident 



in the low mass range, where the formation of large scale 
disks of star- forming gas is highly suppressed (e.g., galax- 
ies g2743 and g85). In contrast, simulations with no winds 
usually form more massive and more concentrated galaxies 
compared to the vzw model, with most of the star formation 
happening in the central regions of galaxies. 



4 THE EFFECTS OF WINDS 
4.1 Time evolution 

Each galaxy is identified not only across simulations with 
different wind models but also back in time at all available 
snapshots, as the central galaxy in the resimulated halo. 
Figure [3] shows examples of the time evolution of galax- 
ies in terms of their stellar mass, gas fraction, SFR, and 
metallicity. Since we are interested in the evolution of galax- 
ies over cosmological time scales, the time evolution of all 
physical quantities has been smoothed over time intervals of 
~ 50 Myr. Major mergers can still be identified as abrupt 
changes in the stellar mass of galaxies (e.g., galaxy g222 be- 
tween redshifts z = 2-3 and g54 at z « 4). These events are 
generally followed by a temporary increase in the SFR of 
galaxies but their overall evolution is dominated by smooth 
accretion, as has been report ed for simulations of high red- 
shift galaxies elsewhere (e.g., iDekel et al.ll2009l : iDave et all 
l201Ch . 

Simulations with no winds clearly result in higher stellar 
masses for all galaxies at all times. Very early on, the SFR of 
galaxies increases very rapidly for simulations with no winds. 
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Figure 3. Time evolution of galaxies g222, g2403, and g54 from z = 8 to z = 2 for simulations including models for momentum-driven 
winds (green solid lines), constant winds (red dashed lines), and no winds (blue dot-dashed lines). From top to bottom: Stellar mass, gas 
fraction, star formation rate, specific star formation rate, and average gas phase metallicity (SFR- weighted). The time evolution of all 
physical quantities has been smoothed over time intervals of ~ 50 Myr. 




R (kpc) 



R (kpc) 



R (kpc) 



Figure 4. Radial profiles of galaxies g222, g2403, and g54 at z = 2 for simulations including models for momentum-driven winds (green 
solid lines), constant winds (red dashed lines), and no winds (blue dot-dashed lines). From top to bottom: Azimuthally averaged face-on 
stellar and gas surface densities, circular velocity (calculated as = GM enc (r)/r for the enclosed mass M enc (r) within radius r), and 
average gas phase metallicity (SFR- weighted). Error bars in the bottom panel indicate the spread of the particle distribution within each 
radial bin for each galaxy. Black vertical dashed lines show the softening length of the simulation at z = 2 (e « 224 pc physical) while 
the green, red, and blue vertical doted lines correspond to the radius enclosing half of the total SFR for the different wind prescriptions. 
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When galactic outflows are included, galaxies regulate them- 
selves resulting in significantly lower SFRs relative to the no 
wind model at z > 3. For momentum-driven winds, recycling 
of the outflowing gas back into galaxies occurs on a time 
scale less than the Hubble time, providing an additional gas 
supply which is not av ailable in simulations with no winds 
(jOppenheimer &; Dave1l201C)h . This additional gas supply re- 
sults in comparatively higher SFRs towards the end of the 
simulation at z = 2. Constant winds, however, are more 
efficient in removing cold gas from galaxies and, therefore, 
result in systematically lower SFRs at all times. Interest- 
ingly, specific SFRs follow a common trend regardless of the 
wind model adopted: galaxies have the highest specific SFRs 
very early on (as high as lOGyr -1 at z = 8) and then de- 
crease monotonically with decreasing redshift. To first order, 
outflows reduce SFRs by a factor 1/(1 +77) relative to halo 
accretion rates (|Dave et al.l l2012h , which leads to specific 
SFRs that remain the same across wind models (SFRs and 
stellar masses are reduced by the same amount). Wind recy- 
cling accretion, however, yields increasing specific SFRs at 
later times, especially in the momentum-driven wind model. 

Gas fractions are systematically lower for the no wind 
simulations compared to either the constant wind or the 
momentum-driven wind models; while galaxies g222, g2403, 
and g54 have gas fractions below 20 % for simulations with 
no winds at z = 2, galactic outflows are able to maintain 
gas fractions a factor of two higher at the end of the simula- 
tion. Gas fractions may be conveniently expressed in terms 
of the depletion time tdep = M^sl SF R and the specifi c 
SFR as / gas = [(1 + (UepsSFR)- 1 )}- 1 J Dave et alj|20l3 ). 
Since specific SFRs are roughly insensitive to wind model 
early on and the depletion ti me is shorter in hig her stellar 
mass galaxies (tdep oc Mr ' 3 : [Dave et aTll2011bh , no wind 
simulations yield lower gas fractions. 

Finally, metallicities are significantly higher for all 
galaxies in the no wind model, owing to a greater conver- 
sion of baryons into stars and a lack of ejection of metals 
into the IGM. These results are qualitatively consistent with 
those obt ained using lower- resolution non-zoom cosmologi- 
cal runs (|Dave et alJl2011aJ lbh. 

4.2 Radial profiles 

The radial profiles of galaxies reveal interesting differences 
among wind models, as shown in Figured for galaxies g222, 
g2403, and g54 at z = 2. Face-on projected stellar and gas 
surface densities, circular velocity, and metallicity as a func- 
tion of radial distance form the centers of galaxies have been 
calculated as average values within logarithmically spaced 
radial bins. The most massive galaxies (g222 and g2403) 
reach resolved stellar and gas surface densities in the range 
S star = 10 4 -10 5 M pc 2 and S gas = 10 3 -10 4 M pc 2 respec- 
tively, with differences up to an order of magnitude among 
wind models. At their centers, stellar surface densities of 
most galaxies reach values well above 10 3 M pc 2 despite 
gravitational forces being softened at scales ~ e « 224 pc 
(Figure El) . 

No wind simulations usually result in higher stellar sur- 
face densities at all radii compared to the vzw and cw wind 
models, and often steeper density gradients at the centers 
of galaxies. The effects of galactic winds are, however, more 
evident in the distribution of the star forming gas: nw galax- 



ies contain most of their gas within the inner kiloparsec, 
reaching gas surface densities which can be factors of a few 
above the cw and vzw models. Galactic winds yield more 
extended disks by removing gas from their centers and hav- 
ing it re-accreted at larger scales. Wind recycling, especially 
in the vzw model, results in significant amounts of gas ex- 
tending up to a few kiloparsecs away from the galaxy center 
(see galaxies g222 and g54 in Figure [4]). As we show in the 
next sections, this result is common to most of our simu- 
lated galaxies. Galaxy g2403 is one of the few galaxies that 
present a very compact gas distribution even in the pres- 
ence of momentum-driven winds, similar to galaxies from 
simulations with no galactic outflows. 

Most of our simulated galaxies are characterized by flat 
rotation curves extending up to a few kiloparsecs. The no 
wind model usually produces more massive and more com- 
pact galaxies and this results in higher circular velocities 
that peak at smaller radii compared to the vzw and cw mod- 
els. Simulations with winds produce more smoothly rising 
rotation curves, with lower circular velocities at small radii 
that tend to approach values similar to the nw model at 
larger scales, where the contribution from the dark matter 
component becomes dominant (see galaxies g222 and g54 
in Figure 31 this is explored in more detail in Section [4.4p . 
Galaxy g2403 is again an atypical case, with very similar 
rotation curves for the nw and vzw wind models. 

Finally, negative radial gas-phase metallicity gradients 
(higher metallicities at the centers of galaxies) are common 
in most simulated galaxies, in general agreement with ob- 
servations (but see lQuevrel et al.ll2012h . Abrupt changes in 
azimuthally averaged metallicities or even positive metallic- 
ity gradients may eventually occur for galaxies undergoing 
mergers (e.g., galaxy g222 at R ~ 3-4 kpc). We find that 
simulations with no winds usually result in higher metallic- 
ities at all radii and also show steeper metallicity gradients. 
Galactic outflows serve to redistribute metal-enriched gas 
from the central regions of galaxies over larger scales, re- 
sulting in less steep metallicity gradients. 

4.3 Connecting individual halos across wind 
models 

In order to understand the global effects of galactic winds it 
is useful to look at their effects on different galaxy properties 
at fixed halo mass, since the latter is expected to be roughly 
insensitive to baryonic processes. In Figure [5j we plot the 
halo baryonic fraction, the central stellar mass to halo mass 
fraction, and the SFR as a function of halo mass for all sim- 
ulated galaxies at z = 2. Here, each vzw galaxy is connected 
to the corresponding nw and cw galaxies to help identify 
both the global effects of outflows on the galaxy popula- 
tion as well as the effects on individual galaxies. Total halo 
masses are very similar for the nw and vzw models at z = 2, 
except occasionally due to the slightly different timing of 
merger events and our adopted definition of dark matter 
halos for each central galaxy (see Section [2]) . Halo masses 
are, however, systematically lower for simulations with con- 
stant winds, for which gas outflows are able to escape more 
easily from the halo potential well, all of which generally 
have escape velocities below the constant wind velocity of 
v w = 680 km s -1 . Note that momentum-driven wind veloc- 
ities are below that of the cw model in all but the most 
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Figure 5. Effects of galactic winds on global properties of z = 2 galaxies and their parent halos as a function of halo mass. Left: Total 
halo baryonic fraction in units of the cosmological baryonic fraction (f^ = 0.165). Middle: Ratio of central stellar mass to halo mass 
in units of f^. Right: Central galaxy star formation rate. Galaxies from simulations including momentum-driven winds (vzw), constant 
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connect galaxies from the vzw simulations with their galaxy counterparts from the cw and nw simulations respectively. 
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massive galaxies, especially a t high redshift when galacti c 
velocity dispersions are lower (|Oppenheimer &; Davel[2008l ). 

Figure [5] shows that halo baryonic fractions are only 
slightly higher for the no wind model (and higher on av- 
erage than the cosmological baryonic fraction = 0.165) 
compared to the vzw model: most of the outflowing gas 



does not escape the halo potential well for the momentum- 
driven wind model. This is despite the fact that the wind 
speed typically is comparable to or exceeding the escape 
velocity (and scales with it as well); but hydrodynamic 
(i.e., ram pressure) slowing is act ually dominant in many 
cases ([Oppenheimer & Dave 2008). In contrast, the high ef- 
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ficiency of constant winds in expelling gas from halos results 
in systematically lower halo baryonic fractions, a significant 
suppression of star formation in the central galaxy, and the 
consequent reduction of stellar masses for all galaxies in our 
sample. 

Figure [5] (middle panel) confirms our earlier expecta- 
tions: for a given halo mass, central galaxy stellar masses 
are systematically higher in the absence of galactic outflows. 
Interestingly, SFRs for the nw and vzw models happen to 
be similar at z — 2 for most simulated galaxies. However, 
as we have seen in Figure [3] for galaxies g222, g2403, and 
g54, momentum-driven winds result in rather different star 
formation histories, with significant suppression of star for- 
mation early on an d enhanced activity due to wind recy- 
cling at later times (jOppenheimer fc Dave 2010). The wind 
recycling channel does not act on the cw simulations for the 
range of halo masses probed here, since these halos generally 
have escape velocities lower than the assumed wind speed. 

Figure [6] illustrates the effects of galactic winds in the 
star formation histories of our simulated galaxies by showing 
the evolution of the cosmic star formation efficiency (cSFE) 
in the redshift range z — 2-6. The cSFE is the ratio of 
the SFR of the centr al galaxy to the to tal halo baryonic 
accretion rate given in Dekel et al. 

Mace ~ 6.6 / b (Mhaio/IO^M©) 1 - 15 (1 + z) 2 - 25 M yr- 1 (1) 

At z — 6, the average cSFE of galaxies for simulations in- 
cluding outflows (cw and vzw models) is ~ 0.15, about a 
factor of 2 lower than simulations with no winds. From z — 6 
to z = 4, the average cSFE increases rapidly up to ~ 0.69 
for nw simulations, and then decreases down to ~ 0.54 at 
z = 2. Simulations with constant winds seem to follow a 
similar tendency, with the average cSFE peaking at z ~ 4- 
5 and decreasing at lower redshifts, but with significantly 
lower values (roughly x3, as expected from the (1 + 77) sup- 
pression of star formation). Momentum-driven winds result, 
however, in an overall increase of cSFEs from z — 6 to z — 2, 
owing to the recycling of gas that was launched into winds 
at earlier times. The vzw model causes an effective delay in 
star formation by ejecting significant amounts of gas from 
small, early galaxies and having it reaccreted by z ~ 2. This 
coincidentally produces average SFRs similar to simulations 
with no winds at z ~ 2. 

4.4 Rotation curves 

Galaxy rotation curves contain substantial information 
about the structure and kinematics of galaxies. Figure 
(top panel) shows the rotation curves of all simulated galax- 
ies at z = 2 (dotted lines), color-coded by wind model. 
Here, we calculate the rotation velocity from the enclosed 
mass at a given radius, which we call v c i rc . Our galaxies are 
characterized by flat rotation curves extending up to several 
kpc scales, with asymptotic circular velocities ranging from 
~ 100 to 260 kins" 1 at z = 2 for the vzw model. In Figure[4l 
we showed for galaxies g222, g2403, and g54 that galactic 
outflows can have a significant impact on the overall shape 
of their rotation curves. 

The thick solid lines in Figure (top panel) show 
stacked rotation curves for each wind model, where circu- 
lar velocities have been averaged over all galaxies as a func- 
tion of radius. Simulations with winds produce galaxies with 
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Figure 7. Top: Circular velocity (i; c irc) as a function of radial dis- 
tance from the center of simulated galaxies at z = 2. We compute 
^circ = ^encW/ r f° r the total enclosed mass M enc (r) within 
distance r. Middle: Azimuthal velocity (v^) of the gas particles 
with respect to the total angular momentum of the galaxy aver- 
aged within logarithmically spaced radial bins (SFR- weighted). 
Bottom: Ratio of the average azimuthal velocity of the gas to the 
circular velocity as a function of radial distance from the center of 
galaxies. Green, red, and blue dotted lines show rotation curves 
of individual galaxies (or the ratio i^/i^irc) for the momentum- 
driven winds, constant winds, and no wind models respectively. 
Thick solid lines show the rotation curve (or v^/v c [ TC ) averaged 
over all galaxies corresponding to each wind model. 

more gradually-rising rotation curves compared to the more 
centrally peaked rotation curves of galaxies with no winds. 
The average circular velocities for the nw and vzw mod- 
els tend to approach similar values at larger radii, where the 
contribution from the dark matter component becomes dom- 
inant. Simulations with constant winds are more efficient in 
removing gas from halos, resulting in systematically lower 
asymptotic circular velocities, but their shape is similar to 
the momentum-driven winds. 

The middle panel of Figure [7] shows a more 
observationally-motivated way of calculating the rotation 
curve. Rather than taking the enclosed mass and assuming 
full rotational support, here we directly calculate the SFR- 
weighted azimuthal velocities of gas particles with respect 
to the total angular momentum of the galaxy, and average 
them within logarithmically spaced radial bins, which we call 
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Figure 8. Build up of the galaxy rotation curves as a function 
of redshift from z = 6 (thin grey line) to z = 2 (thick black line) 
for simulations including momentum-driven winds. Dashed lines 
show the circular velocity (i^irc) averaged over all galaxies (as 
in Figure [3 top panel) while solid lines show the actual gas az- 
imuthal velocity (v^) averaged within each radial bin (and again 
averaged over all vzw galaxies). 



V(/>. This is analogous to an Ha rotation curve. The rotation 
curves derived in this way show a significantly more gradual 
increase in circular velocity with radius in the inner parts, 
because of the increased dispersion support in the central 
regions. This is particularly clear for simulations with no 
winds, which produce more compact galaxies and centrally 
peaked rotation curves based on the mass profiles. 

The bottom panel of Figure [7]shows the ratio of v^/vdrc- 
We see that the gas azimuthal velocity is, on average, lower 
than the inferred circular velocities. The gas reaches rota- 
tion velocities comparable to the circular velocity at ^ 3 kpc 
away from the centers of galaxies. In the inner regions, the 
contribution of random motions to the dynamical equilib- 
rium of galaxies becomes comparable to the support from 
ordered rotation. The discrepancy between and v c \rc at 
scales < 3 kpc occurs for all simulated galaxies and wind 
models, and it is well resolved in our simulations. Using 
to infer the mass profile of galaxies could, therefo re, lead to 
significant misestimations ( Valenzuela et al"1l2007l ). 

Figure [8]shows the build-up of rotation curves with time 
by comparing stacked rotation curves from z = 6 to z = 2 for 
our simulations with momentum-driven winds. Solid lines 
show Vc/) (for the gas), and dashed lines show v c irc- We find 
that the average asymptotic circular velocities increase from 
~ 90 to 170 km s -1 with decreasing redshift, as expected be- 
cause our galaxies are becoming more massive. For v c irc, the 
peak of the (average) rotation curve is always at ~ 1 — 2 kpc, 
perhaps moving inwards to lower redshift s. In contrast, the 
peak of is at 3 — 4 kpc, and does not vary much with 
redshift. Comparing to v c irc, we see that gas azimuthal 
velocities are, on average, lower than the inferred circular 
velocities at all redshifts. The ratio v^/vdrc evaluated at 
R ~ 3 kpc increases from ~ 0.48 at z = 6 to ~ 0.86 at 
z = 2, indicating that galaxies become, on average, progres- 
sively more rotationally supported with time (and therefore 
with increasing mass) in the redshift range z = 6 — >• 2. 



5 COMPARISON TO OBSERVATIONS 

The high resolution afforded by Hubble probing both young 
stars and older stellar populations, together with adap- 
tive optics-enhanced spectral studies from the ground, have 
opened up new windows into high redshift galaxy assembly. 
Here we examine both qualitatively and quantitatively how 
our simulated galaxies compare with the latest observations 
of z ~ 2 galaxies. 



5.1 Qualitative galaxy morphologies 

We begin by considering in detail the morphology of the two 
galaxies that were specifically chosen for resimulation, g54 
and g222, and therefore lie near the center of the zoomed 
region. Recall that g54 was chosen to be somewhat more 
isolated and have a fairly quiet merger history, while g222 
lived in a denser region with a more violent merger history. 

Figure [9] shows two-dimensional projected views of 
galaxy g54 at four different inclination angles. The top three 
rows show the stellar surface density, gas surface density, and 
SFR density distributions resulting from our fiducial simu- 
lation including momentum-driven winds. The three bottom 
rows show mock Ha line intensity, line-of-sight velocity, and 
velocity dispersion maps of galaxy g54 that mimic the res- 
olution of current near-infrared integral field spectroscopic 
observations with SINFONI. Simulated SFRs have been con- 
verted to Ha lumi nosity using Lhol [erg s -1 ] = 2.3 x 10^_SFR 
[M Q y r -1 ] (from lKennicuttlfl998bl corrected for a Cha brierl 
(2003) IMF). Then, mock Ha line intensity maps have been 
obtained by (i) placing the simulated galaxy at z = 2, con- 
verting its physical size to its apparent angular size, (ii) 
convolving the obtained Ha flux map with a 0.1 7" beam, 
and (Hi) matching the pixel size to those of typical obser- 
vations (0.05"/pix). Mock Ha flux- weighted line-of-sight ve- 
locity and velocity dispersion maps have been obtained by 
degrading the spatial resolution in a similar way. 

With stellar mass ~ 2.5 x 10 10 M©, gas fraction / gas ~ 
0.48, and SFR ~ 13.6 M© yr -1 , galaxy g54 has a prominent 
two-arm spiral structure extending up to ~ 10 kpc scales 
that could be observable even at high inclination angles 
(i ^ 60°). Star forming gas and stars generally trace the 
same structure, though the stars show a prominent bulge 
that is absent in the gaseous component. Mock Ha line in- 
tensity maps show a somewhat compact central nucleus with 
radius ~ 1 kpc and a factor ~ 2 increased brightness relative 
to the spiral arms (depending on inclination angle). Obser- 
vations using star formation tracers and rest-frame optical 
observations should, thus, infer similar spiral morphologies. 
The velocity maps illustrate the effects of the inclination an- 
gle on the observed kinematic properties of galaxies. A very 
smooth velocity gradient along the morphological major axis 
can be clearly identified for this galaxy at a wide range of 
inclination angles (i ^ 30°). Interestingly, galaxy g54 shows 
a characteristic "spider diagram" pattern in the velocity iso- 
contours at intermediate inclination angles (i « 30° — 60°), 
as expected for inclined rotating disks. For nearly edge-on 
observations, line-of-sight projected peak-to-peak velocities 
reach values > 400 km s -1 , consistent with the circular ve- 
locity Vcirc ~ 230kms _1 calculated at the outer edge of the 
disk (R ~ 10 kpc). With a relatively flat velocity dispersion 
map, this galaxy would be identified as a large, thick, rota- 
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Figure 9. Morphological and kinematic properties of galaxy g54 at 2 = 2 for our fiducial simulation including momentum-driven winds. 
Upper panels show (from top to bottom) the projected stellar surface density, gas surface density, and SFR surface density color coded 
by (SFR-weighted) line of sight velocity. Lower panels show (from top to bottom) mock Ha line intensity, line of sight velocity, and one- 
dimensional projected velocity dispersion maps that mimic integral field unit observations with SINFONI. Two-dimensional projected 
images and mock maps are shown at different inclination angles, from direct face-on (left) to edge-on (right). The region shown is 20kpc 
(physical) across. 
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Figure 10. Same as Figure [9] for galaxy g222. 



tion dominated disk. These properties are in general quite 
similar to, though somewhat scaled dow n from BX442, th e 
z — 2.18 grand design spiral observed by lLaw et al.l (|2012f ). 
BX442 has a measured stellar mass ~ 6 x 1O 1O M , Ha- 
derived SFR ~ 45M0yr _1 , inclination-corrected circular 
velocity v c irc ~ 234kms _1 , and very similar size and mor- 
phology compared to our simulated galaxy g54. 

Figure [10] shows similar two-dimensional maps for 



galaxy g222, a marginally higher mass galaxy that has un- 
dergone more frequent interactions and mergers. Indeed, 
there is an infalling galaxy along the lower spiral arm, which 
is prominent in the stellar distribution but much less evident 
in the gas, possibly because its gas has been stripped during 
infall. The stellar surface density distribution shows a high 
concentration of stars at the center of g222 and possibly a 
weak, extended spiral structure. The gas and SFR distri- 
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butions reveal high levels of turbulence on this galaxy. The 
morphology of g222 appears highly disturbed in Ha emission 
and it does not seem to trace the smooth stellar component 
as in the case of g54. The velocity maps show signs of or- 
dered rotation in the underlying large scale gas distribution 
at high inclination angles (i ^ 60°), but the inferred velocity 
gradients are significantly disturbed. The velocity dispersion 
map reveals an irregular gas clumpy structure with signifi- 
cant turbulent motions. The Ha shows up to several clumps, 
depending on the viewing angles, but these clumps are not 
apparent in the stellar distribution, which suggests these are 
shor t-lived gaseous cl umps as found in similar simulations 
by iGenel et all (l2012h . and as inferred in observations by 
IWuvts et al. (|2012h ; we will examine clump properties in 
more detail in future work. In general, this simulated galaxy 
is b roadly similar to clumpy disk galax i es in the SINS sam- 
ple (|Forster Schreiber et alJl2009U201ll ; lGenzel et al.ll201lh . 

These two examples illustrate the wide diversity in mor- 
phologies at z — 2 predicted for disk galaxies in these simu- 
lations, as well as its dependence on tracer (gas, stars, Ha). 
Disk structures are usually present, but can range from qui- 
escent "grand-design" spirals to turbulent and clumpy disks. 
Next, we quantitatively compare basic structural properties 
of our full sample of simulated galaxies to recent observa- 
tional surveys of z ~ 2 disks. 

5.2 Comparison to SINS and SHiZELS 

Our galaxies appear to be qualitatively similar to those ob- 
served at z rsj 2. We now conduct a more quantitative com- 
parison to the structural and dynamical properties of z ~ 2 
galaxies inferred from near-infrared integral field spectro- 
scopic observations obtained with SINFONI at the VLT. In 
particular, we fo cus on the "Ha sample" of th e SINS sur- 
vey presented in iForster Schreiber et al.l (|2009h , consisting 
of 62 rest-frame UV/optically selected star forming galaxies 
at z & 1.5-2.5, and the SHiZELS survey ([Swinbank et alj 
l2012h . consisting of 9 Ha-selected galaxies at z w 0.84-2.23 
drawn from the H iZELS near-infrared narrow-band survey 
(jGeach et alJ l2008). In Figure [HQ we show some key struc- 
tural and dynamical quantities versus stellar mass, for our 
simulated galaxies color-coded by wind model, as well as for 
the SINS (black crosses) and SHiZELS (brown diamonds) 
samples. 

The upper left panel of Figure [11] shows the specific 
SFR as a function of stellar mass for our sample of galax- 
ies at z = 2. Our simulations naturally predict the exis- 
tence of a correlation between galactic SFRs and stellar 
masses, M* oc SFR, fairly insensitive to stellar feedback 
models, due to the dominance of smooth and steady cold 
accret ion (e.g., iFinlator et aT]|2006l ; iDave 2008; D ave et all 
l2011ah . As expected from models, a tight M*-SFR relation 
has been observed out to z ~ 2, with a roughly constant 
slope close to unity (e.g ., iDaddi et al.l 120071 ; lElbaz et al.l 
2007; iNoeske et~aH 120071 ). However, the observed normal- 
ization of this relations increases more rapidly with red- 
shift from z = —> 2 that current cosmological simula- 
tions pre dict, resulting in too low pred icted specific SFRs 
at z ~ 2 (IDaddi et al.ll2007l : lDavell2008l ). Our zoom simula- 
tions follow this trend - Fig ure | 11 | shows that t hey lie below 
the M*-SFR correlation of IDaddi et al.1 (|2007T ) by a factor 
x2 — 3 in specific SFRs (for the vzw model), and by some- 



what higher factors for the other wind models. The discrep- 
ancy is slightly worse when compared to the SINS galaxies, 
although the Ha selection could pla y some role in preferen- 
tially picking out high-SFR objects (|Forster Schreiber et al.l 
2009). On the other hand, the difference between our sim- 
ulations and the SHiZELS galaxies appears to be not as 
large. This discrepancy between models and data has been 
much debated; for instance recent Herschel results suggest 
that Ha-inferred SFRs may slightly overestimate the true 
(i.e., far-infrared bolom etric) SFR by some non-trivial fac- 
tor (|Nordon et al.ll2010l ) . Other ideas for explaining this dis- 
crepa ncy invoke variations in t he stellar initial mass func- 
tion (Naravanan fe Pavel [2012). We will not pursue this 
here, except to note that our zoom simulations do not alle- 
viate this discrepancy known from our larger cosmological 
runs. 

Simulations with momentum-driven winds result in 
higher specific SFRs compared to simulations with no winds 
at z = 2, more in agreement with observations. In the vzw 
model, the mass loading factor scales as 77 ~ 1/a and out- 
flows tend to suppress early star formation while providing 
high gas fractions to maintain comparatively high SFRs at 
later (z ~ 2) times, as seen in Figures [3] and [6] This suggests 
that even stronger outflows might be needed i n simulation s 
in order to match the z ~ 2 M*-SFR relation (jPavell2008h . 
This is not trivial, however, since outflows must not be too 
strong lest th ey fail to produce enough early m etals to en- 
rich the IGM ( Qppenheimer fe Pavell2006lJ2008l) and enoug h 
photons to reionize the Universe (e.g. jFinlator et al1l2012h . 

The middle left panel of Figure [TT1 shows the z = 2 size- 
mass relation for all simulated galaxies. The effective radius 
-Ri/2, s fr is defined to enclose half of the total SFR and it is 
calculated as a two-dimensional projected radius averaged 
over 100 random viewing angles. Here we see that there is a 
clear separation of galaxies in the M*— -Ri/2 jS fr plane for the 
different wind models. For a given stellar mass, simulations 
with no winds produce galaxies with a very compact distri- 
bution of gas, with most of the star formation happening 
within < 1 kpc from their centers. Simulations with con- 
stant winds populate the low mass and small size region of 
the plot, and only the more massive galaxies in the sample 
are able to maintain gas disks with sizes i?i/2, s fr > 1 kpc. 
Momentum-driven winds produce more extended galaxies 
with sizes of several kpc, along with a minority of more 
compact galaxies. 

Outflows affect sizes by preferentially ejecting the star- 
forming gas from t he central regions and having it re- accrete 
over larger scales ([Brook et al.l l2012), resulting in large, ex- 
tended star-forming disks that are more in agreement with 
SINS and SHiZELS data. Nonetheless, it is interesting that 
compact galaxies can also occur in the vzw case; galaxy 
g2403 is one such compact galaxy (see Figure 2}. From our 
visualization^], it appears that these more compact systems 
gener ally arise during or sh ortly after a significant merger 
event (Bou rnaud et al ]|201l|), but we will quantify this more 
rigorously in the future. 

The diversity of sizes is also reflected in the half-mass 
radii of the stars R1/2,* (Figure [TT1 middle right). Generally, 
the sizes are comparable to the half-SFR radius, but in some 



1 See http://www.physics.arizona.edu/~angles/movies/ 
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Figure 11. Structural and dynamical quantities calculated from our sample of simulated galaxies at z = 2, for our momentum-driven 
wind (vzw; green stars), constant wind (cw; red triangles), and no wind (nw; b lue squares) simulations. Blac k crosses with grey error bars 
in each panel represent the z ~ 2 star-forming galaxies of the SINS survey (Forster Schreiber et al 1 l2009h . while the brown diamonds 
correspond to a sample of galaxies from the SHiZELS survey (|Swinbank et al.ll2012l ). Toy left: Specific star forma tion rate (sSFR) as a 
function of stellar mass (M*). The shaded beige band shows the z = 2 M*-sSFR correlation of Daddi et al. ( 2007) with a scatter (semi- 
interquartile range) of 0.16 dex in sSFRs. Top right: Stellar bulge mass fraction as a function of M*. Middle left: Radius enclosing half of 
the total star formation rate as a function of M*. Middle right: Radius enclosing half of the total stellar mass as a function of M*. The 
dashe d line shows the nominal threshold stellar surface density M*/i?J a 5 lf = 10 10-3 Mq kpc — 1-5 for compact galaxies from lBarro et al] 
(|2012h . Bottom left: Velocity dispersion a of the gas component of simulated galaxies as a function of M*. Bottom right: Disk peak 
rotational velocity divided by a of the star forming gas, as a function of M* for all simulated galaxies. 
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cases R1/2,* can be < 1 kpc. This may have interesting con- 
sequences for the progenitors of so-called compact ellipticals, 

1. e., early- type galaxies at z ~ 2 that show very high stel- 
lar densities. We show as the dashed line a nominal thresh- 
old stellar surface density of MjR]f = 10* 03 M kpc -15 
for compact galaxies from iBarro et all (|2012l , see their Fig- 
ure 1). The radius used for these observations is actually 
the iif-band half-light radius, but this should be fairly com- 
parable to the stellar half-mass radius. Even with winds, 
a small fraction of our galaxies lie above this surface den- 
sity threshold (i.e., below the line). This suggests that our 
simulations do, with some frequency, produce galaxies that 
have sufficient stellar densities to be the progenitors of com- 
pact ellipticals. These simulations have no mechanism for 
quenching star formation, and hence our compact galaxies 
are still star-forming; it is an additional constraint on mod- 
els to form compact passive systems that whatever feedback 
mechanism is responsible for quenching, it operates when 
the galaxy is in a dense (likely post-merger) state. 

Galactic outflows seem to affect the bulge mass frac- 
tion of galaxies, as shown in Figure 111! upper right panel. 
To estimate the stellar bulge mass fraction, we perform a 
simple bulge-disk kinematic decomposition for the stellar 
content of all simulated galaxies. We calculate the azimuthal 
velocity of each star particle with respect to the direc- 
tion of the total angular momentum of the galaxy as for 
our rotation curves, and estimate the bulge mass as double 
the mass of particles moving with < 0. Though this is 
not analogous to observational determinations of the bulge 
mass, it accurately characterises the mass in the spheroidal 
component of our simulated systems. With winds, the bulge 
fraction decreases with increasing stellar mass, suggesting 
that galaxies start out small and dispersion-dominated, 
and move towards being more ordered disks. The opposite 
trend is seen for our no wind simulations, tracking more 
the canonical behavior that galaxies start out as disks and 
then merg e together to form mor e dispersion-dominated sys- 
tems fe.g. JWhite &; Frenklll99lh . SINS observations suggest 
that smaller galaxies tend to be more dispersion-dominated, 
qualitatively favoring our wind simulations, although a more 
careful comparison that accurately mimics how bulge-to- 
disk ratios are measured in data are needed for a more defini- 
tive result. 

The lower left panel of Figure [11] shows the gas ve- 
locity dispersion, a, as a function of stellar mass at z = 

2. Here a is calculated as the spatially integrated (SFR- 
weighted) one-dimensional velocity dispersion calculated 
within i?i/2 >s fr an d averaged over 100 random orienta- 
tions. Despite being extended and rotationally supported 
disks, most of our galaxies are characterized by high ve- 
locity dispersions (> 30kms -1 ). This suggests that tur- 
bulent motions are significant even in rotationally sup- 
ported z ~ 2 galaxies, as inferred from observations 

Law et al.ll2009l:lForster Schreiber et al.ll2009l : IWright et all 



2009; ISwinbank et al.|[2Q12h and reported in previou s simu- 
lations (e.g. JCeverino et al.ll201Ct iGenel et al.ll2012i ). For a 
given stellar mass, simulations including outflows result in 
galaxies with higher gas velocity dispersion, in better agree- 
ment with observations. Note that since outflowing gas is hy- 
drodynamically decoupled as it is ejected from the ISM, the 
outflows themselves are not directly injecting turbulence. In- 
stead, this increased turbulence relative to simulations with- 



out winds owes to the higher gas fractions and the recycling 
of the outflowing gas back into the galaxies. 

The lower right panel of Figure ITTl shows the disk rota- 
tion velocity, Vd divided by velocity dispersion a, as a func- 
tion of stellar mass. Here Vd is taken as the peak azimuthal 
velocity (v^) from the gas rotation curves obtained in Sec- 
tion [431 The momentum-driven wind model yields a typical 
Vd/cr ~ 3, which is in broad agreement with the turbulent 
high-z disks seen in SINS and SHiZELS. Comparisons with 
data are not straightforward since rotation is computed from 
a variety of means in the data, and we have not tried to 
mimic this in detail; hence we have not plotted data here. 
Constant and no wind models have slightly more rotational 
support on average, though still within the range of the data. 
There is little trend with mass, much less than for the bulge 
fraction, so even though more massive galaxies are diskier, 
they are not more rotationally-supported. 

Overall, the morphologies and dynamical properties of 
simulated disks at z — 2 are in fair agreement with obser- 
vations from the latest integral field unit studies of high-z 
star- forming galaxies, particularly in the case of momentum- 
driven winds. This occurs despite some overly simple as- 
sumptions in the modeling, such as decoupling of wind ma- 
terial escaping the disk, and the usual concerns about the 
ability of SPH to s uppress viscosity a nd resolve dynami- 
cal instabilities fe.g. jAgertz et al]|2007h . This suggests that 
either these neglected/poorly modeled processes are not im- 
portant for establishing the structural properties of disks, 
or else these agreements with data are merely coincidental. 
Nonetheless, these simulations appear to provide at least a 
plausible model for the formation of disks during this cosmic 
epoch. 

5.3 Evolution of Physical Properties 

Using our plausible simulated population of z = 2 disk 
galaxies, we now examine the evolution of physical prop- 
erties from z — 8 — » 2. We focus here on the momentum- 
driven wind model, since from a variety of wider constraints 
it is our favoured mod el of the three presented here (see, 
e.g., iDave et alJkoilal lbh. However, many of the results are 
broadly similar for our other wind models. We show the 
evolutionary tracks for our 8 galaxies in various colours in 
Figure [121 Since we are interested on the evolution of galax- 
ies on cosmological time scales, all physical quantities have 
been averaged over time intervals of ~ 150 Myr. All tracks 
go from left to right (i.e., lower to higher mass), with indi- 
vidual unit redshifts indicated by the points along the tracks 
starting at z = 8. The same data as shown in Figure ITTl is 
reproduced here, but note that this is for observed z ~ 2 
galaxies and hence it is shown here for reference. 

The top left panel of Figure [12] shows the evolution- 
ary tracks of simulated galaxies in the main sequence (M*- 
sSFR) plane. All galaxies show a similar evolution of de- 
creasing specific SFR as they increase their stellar mass, 
with a steeper slope relative to the observed z = 2 star 
forming main sequence. To first order, this is consistent 
with a roughly linear SFR to M* relation, with the over- 
all normalization decreasing with redshift in a way simi- 
lar to the total mass-accretion rate on to haloes, ~ (1 + 
z) 2 25 , as predicted from cold accretion- driven galaxy growth 
(|Dekel et al.1 12009). Momentum-driven winds make galax- 
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Figure 12. Evolution of structural and dynamical quantities calculated from our sample of simulated galaxies from z = 8 — > 2. Different 
color tracks show each of our eight momentum-driven wind run galaxies, with filled circles indicated at integer redshifts (stars show 
the location of galaxies at the end of the simulation at z = 2). Other lines and symbols are as in Figure [TT] Top left: Specific star 
formation rate (sSFR) as a function of stellar mass (M*). Top right: Stellar bulge mass fraction as a function of M*. Middle left: Radius 
enclosing half of the total star formation rate as a function of M*. Middle right: Radius enclosing half of the total stellar mass as a 
function of M*. Bottom left: Velocity dispersion a of the gas component of simulated galaxies as a function of M*. Bottom right: Disk 
peak rotational velocity divided by a of the star forming gas, as a function of M* for all simulated galaxies. The time evolution of all 
physical quantities has been averaged over time intervals of ~ 150 Myr. 



ies deviate from this characteristic cosmological inflow rate 
by suppressing star formation at early epochs and recycling 
that material back into galaxies at later times, yielding a 
shallower evolution of the specific SFR with time. 

At z — 2, galaxy sizes scale with stellar mass for sim- 



ulations with winds, in agreement with observations (Fig- 
ure [TT]) . Middle panels in Figure [12] show that galaxies tend 
to increase in size with time, as expected, but their evolu- 
tionary tracks exhibit significant variation between galax- 
ies, resulting in a large scatter in the size-mass diagram 
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at any given redshift. Major mergers, despite representing 
only a fraction of the total mass growt h in galaxies (e.g., 
iKeres et" al. 2005; Rodighiero et alJuOllI ). have a significant 
impact on galaxy sizes, resulting in large variations of the 
half-SFR radius i?i/2, s fr- Major mergers are also evident in 
the evolution of the half-mass radius, R1/2,*, and may drive 
galaxies towards the region of the size-mass diagram popu- 
lated by compact ellipticals. Interestingly, our most compact 
galaxy at z = 2 (g2403) was also in the compact regime at 
z ~ 4. At the earliest epochs, the radius of galaxies does 
not change much, and may in fact become smaller with 
time. This partly reflects our star formation criterion of 
tih > 0.13 cm -3 that results in star formation occuring over 
a wide area when the universe is very dense; it is therefore 
unlikely to be a robust prediction. 

The upper right panel of Figure [12] shows the stellar 
bulge fraction evolution. Early on, most galaxies are bulge- 
dominated, as a result of high merger rates and large gas 
reservoirs that keep the mass distribution disordered. As 
time proceeds, the combined effects of smooth gas accreting 
into galaxies with higher specific angular momentum and 
galactic outflows removing preferentially low angular mo- 
mentum gas from their centers, cause galaxies to increase 
their sizes and reduce their stellar bulge fraction s on average 
(jGovernato et alJl2009l . l2010l : iBrook et alJfepllh . It is inter- 
esting that this process only kicks in around z ~ 4 to start 
producing disk-dominated systems, at least for the range of 
galaxy masses considered here. Our simulations, thus, pre- 
dict that z ^ 4 was the beginning of the epoch of disk forma- 
tion for massive galaxies. By z ~ 2, this results in a trend of 
decreasing bulge fraction with increasing stellar mass in this 
wind model. We note that this trend is opposite to simple 
expectati ons from classic hiera rchical galaxy formation mod- 
els fe.g. JWhite & Frenk|[l99lh . in which disks form first and 
then merge later to give rise to more dispersion-dominated 
systems. 

The bottom left panel of Figure [12] shows how simulated 
galaxies evolve with redshift in the M*-cr plane. Galaxies 
tend to evolve along the correlation of a and M* observed 
at z ~ 2, extrapolated to lower masses. Disk peak rotation 
velocities follow a similar trend with stellar mass, giving rise 
to a nearly constant ratio Vd/cr « 3 for all simulated galax- 
ies independent of redshift, as shown in Figure [12] bottom 
right panel. Hence while the predicted Vd/cr values are simi- 
lar to that observed at z = 2, it remains to be seen whether 
these same simulations can produce very thin disks as seen 
at z = 0. In order to do so, something must alter the cur- 
rent evolution of Vd/cr, perhaps as a result of the dropping 
accretion rate or else the dropping outflow rate. 

Overall, the redshift evolution of these properties shows 
a consistent buildup of size and velocity dispersion with time 
(and mass), although mergers can particularly impact the in- 
ferred sizes substantially over short time periods. The large 
scatter in observed sizes may, therefore, reflect the short- 
term merger history of galaxies even more reliably than in- 
stantaneous SFRs or bulge-to-disk ratios. Galaxies evolve 
from being bulge-dominated when small to disk-dominated 
when larger, with disks becoming prominent only at z < 4. 
Despite this, galaxies are always rotationally-supported even 
at early times when bulge-dominated, suggesting that the 
present-day association between small bulge-to-disk ratio 



and large rotational support does not necessarily apply to 
high-redshift galaxies. 



6 RESOLUTION CONVERGENCE 

A key computational issue we face in this study is that 
the observations we are comparing against are approach- 
ing the numerical resolution limit of our simulations. Note 
that our nominal ~ 220 pc resolution at z — 2 corresponds 
to the equivalent Plummer force softening length; the scale 
at which we are computing exact gravitational forces is, 
thus, 2.8 times this length, or ~ 600 pc. Furthermore, SPH 
is known to have issues with excessive viscous transport 
of angul ar momentum , and indeed the moving mesh code 
Arepo (ISpringelll2010h is seen to produc e larger-scale disks 
than Gadget ([Vogelsberger etlrt1l2012h from identical ini- 
tial conditions. It is, therefore, important to carry out a basic 
test of resolution convergence of our key results. 

We run simulations with a factor of 2 lower spatial reso- 
lution and a factor of 8 lower mass resolution (i.e., equivalent 
to 512 3 resolution in a [24/i _1 Mpc] 3 box) from otherwise 
identical initial conditions, and identify the galaxies that 
correspond to our high resolution galaxy sample. Figure [13] 
shows the effects of numerical resolution on key structural 
and dynamical quantities of galaxies versus stellar mass for 
our momentum-driven winds and no wind simulations at 
z = 2 (as in Figure [TT}. Here, solid lines (green and blue for 
vzw and nw respectively) connect individual galaxies corre- 
sponding to the higher and lower resolution simulations to 
help identify any systematic trends. 

We note that a galaxy- by- galaxy comparison represents 
a strong numerical convergence test given that even small 
deviations of orbital parameters in galaxy mergers may re- 
sult in rather different structural and kinematic properties 
of remnant galaxies at a given time. Global properties of 
galaxies are thus expected to exhibit better numerical con- 
vergence. Indeed, Figure [13] shows that total stellar masses 
and specific SFRs of galaxies are both well converged for the 
vzw and nw simulations. 

The size of galaxies is somewhat more sensitive to reso- 
lution, as shown in Figure ITTJl middle panels. For individual 
galaxies, we find that the half-SFR (-Ri/2, s fr) an d half-mass 
(^1/2,*) radii may differ by up to a factor of ^ 2 for the 
effective 1024 3 and 512 3 simulations at z = 2. Such sensitiv- 
ity might be expected given that radii track recent merger 
activity, which can vary substantially owing to the chaotic 
nature of orbits within hierarchically growing halos. Despite 
this, we find no systematic trend in the whole sample with 
resolution or wind model - our lower resolution simulations 
confirm that outflows generally produce more extended star 
forming disk galaxies but that may also result in compact 
galaxies (similar to simulations with no winds) with some 
frequency. 

Stellar bulge mass fractions are also slightly sensitive to 
resolution (Figure [13] upper right), but more so for models 
without winds. In the no wind case, bulge fractions are gen- 
erally (though not always) higher at high resolution. Such a 
trend is not as apparent in the vzw case, though in two cases 
the bulge fractions are substantially different. Overall, our 
general result that bulge fractions are higher in simulations 
with no winds remains unaffected by resolution. 
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Figure 13. Resolution convergence of physical quantities at z = 2 for our momentum-driven wind (vzw; green) and no wind (nw; blue) 
simulations, comparing our effective 1024 3 runs (green stars and blue squares for vzw and nw respectively) with effective 512 3 runs 
(green and blue filled circles for vzw and nw respectively) of the same initial conditions. Green and blue solid lines connect individual 
galaxies corresponding to the high resolution and lower resolution simulations for the vzw and nw simulations respectively. Other lines 
and symbols are as in Figure [TT] Top left: Specific star formation rate (sSFR) as a function of stellar mass (M*). Top right: Stellar bulge 
mass fraction as a function of M*. Middle left: Radius enclosing half of the total star formation rate as a function of M*. Middle right: 
Radius enclosing half of the total stellar mass as a function of M* . Bottom left: Velocity dispersion a of the gas component of simulated 
galaxies as a function of M*. Bottom right: Disk peak rotational velocity divided by a of the star forming gas, as a function of M* 
for all simulated galaxies. 
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We find a systematic trend for galaxy velocity disper- 
sion (a) to decrease with resolution (Figure [T3l lower left), 
though again this is significantly more evident for simula- 
tions without winds. This suggests that our effective 1024 3 
simulations may not have reached numerical convergence 
and that higher resolution simulations with outflows could 
result in slightly higher a values for a given stellar mass, 
more in agreement with observations. Despite this, gas recy- 
cling for simulations with outflows seems to result in increas- 
ing a values and, therefore, turbulence in galaxies regard- 
less of numerical resolution. Figure [131 lower right, shows 
that the ratio Vd/cr is reasonably well converged for the vzw 
model but increases systematically for nw simulations with 
lower resolution due to the overall decrease in a. 

In short, our results are not perfectly resolution con- 
verged, and higher resolution generally results in higher 
bulge fractions, smaller sizes, and lower velocity dispersions. 
However, the resolution convergence exhibited by the wind 
simulations is considerably better than that for no winds. 
Interestingly, this trend is also seen for globa l mass and star 
formation rate functions ([Dave et al.ll2011ah . By increasing 
sizes and having a physically-motivated driver of wind recy- 
cling that sets the velocity dispersion, winds actually seem 
to help resolution convergence somewhat. Nonetheless, this 
convergence experiment only spanned a factor of 2 in spa- 
tial scale and 8 in mass, and hence simulations with greater 
dynamic range will be needed to fully assess the numerical 
robustness of these results. 



7 SUMMARY AND CONCLUSIONS 

Powerful galactic outflows are ubiquitous in high redshift 
galaxies and likely play a central role in the evolution of 
galaxies and the IGM. Here, we have presented high resolu- 
tion cosmological zoom simulations that follow the evolution 
of a sample of eight central galaxies down to z = 2, focus- 
ing on the impact of strong outflows on their morphologies, 
kinematics, and star formation properties. Our main results 
are generally resolution-converged and can be summarized 
as follows: 

• Despite the limited sample of eight galaxies presented 
here, our simulated systems span a wide range of morpho- 
logical characteristics at z = 2. Disk structures are prevalent 
but can range from very compact gas and stellar distribu- 
tions, to extended quiescent "grand-design" spirals, to tur- 
bulent and clumpy disks. Inferred morphologies can depend 
on the observed tracer. 

• Simulations with no winds produce rapidly rising SFRs 
that result in higher stellar masses, higher metallicities, and 
lower gas fractions for all galaxies at all times relative to 
simulations with galactic outflows. Momentum-driven winds 
cause an effective delay in star formation by ejecting signif- 
icant amounts of gas from small, early galaxies and hav- 
ing it reaccreted at later times, resulting in higher gas frac- 
tions and star formation histories more in agreement with 
observations. All wind models fail, however, in reproduc- 
ing the normalization of the observed z — 2 M*-SFR re- 
lation, though the late-time recycling by the momentum- 
driven wind model comes the closest; this suggests that even 
stronger outflows at early times and/or small masses may be 
required. 



• Galactic outflows affect the amount and distribution 
of metals in galaxies by regulating star formation, eject- 
ing metals into the surrounding gas preferentially from their 
centers, and the recycling of enriched gas back into galax- 
ies over larger scales. This results in lower metallicities and 
less steep metallicity gradients relative to simulations with 
no winds. The resulting central metallicities are somewhat 
super-solar. Examples of inverted metallicity gradients are 
uncommon among our galaxy sample. 

• Winds have a significant impact on the structural prop- 
erties of simulated galaxies. No wind simulations generally 
produce more compact galaxies with higher stellar surface 
densities, higher stellar mass bulge fractions, and most of 
the star-forming gas concentrated within the inner kilopar- 
sec. Galactic winds usually yield more extended disks and 
tend to reduce bulge fractions by preferentially removing 
low angular momentum gas from their centers. Nonethe- 
less, simulations with winds may produce, in some cases, 
galaxies with stellar surface densities above the threshold 
for compact ellipticals. Sizes in general are quite sensitive 
to the merger history, moreso than star formation or bulge 
fraction, and hence this may be the best way to assess the 
hierarchical buildup of galaxies. 

• Simulations with winds produce galaxies with more 
gradually-rising rotation curves compared to the more cen- 
trally peaked rotation curves of galaxies without winds. 
When calculated from the gas rotation velocities rather than 
the enclosed mass, rotation curves are more smoothly rising 
because of the increased dispersion support in the central 
regions of galaxies. 

• Peak rotation velocities and velocity dispersions scale 
with stellar mass for all wind models, in a manner broadly 
consistent with observations. The inferred ratios Vd/cr are 
consistent with rotationally-supported turbulent disks at 
z = 2. Gas recycling and the high gas fractions of galax- 
ies from simulations with outflows yield higher gas velocity 
dispersions and more turbulent disks compared to no-wind 
simulations. Early small galaxies have both high bulge frac- 
tion and are rotation-dominated, counter to typical trends 
among local galaxies, and suggests that the standard intu- 
ition from today's Hubble sequence may not apply to high- z 
galaxies. 

Our simulations complement previous studies by ana- 
lyzing the effects of large scale outflows on the internal struc- 
ture and evolution of individual z = 2 galaxies. We empl oy 
the same outflow mechanisms used by I Dave et alj (|201 lal lbh 
in non-zoom cosmological simulations, including an obser- 
vationally constrained prescription for momentum-driven 
winds, and with no further tuning of model parameters. 
Encouragingly, simulations with momentum-driven winds, 
which are favoured by a wide ra nge of observations a nd re- 
cent idealized galaxy simulations ([Hopkins et al.ll2012h , yield 
similar trends on the global properties of galaxies when ap- 
plied to simulations with a factor xl50 increased in mass 
resolution, and produce galaxies at z = 2 with structural 
and kinematic properties in broad agreement with observa- 
tions. This provides a new and non-trivial test of hierarchical 
galaxy formation models. 

It is interesting that even among our limited sam- 
ple, we produce a diversity of morphologies, from grand 
design spirals to clumpy turbulent disks to potential 
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progenitors of compact ellipticals. This diversity is en- 
hanced by feedback, and may in fact be governed by it. 
While there is an overall consensus on the importance 
of feedback associated with star formation in determin- 
ing the morphologies of galaxies and their evolution over 
cosmic time, major uncertainties remain on the nature 
and relative significance of different feedback processes. 
The specific outcomes of particular feedback models of- 
ten depend on numerical resolution, the treatment of ISM 
physics, the star fo rmation prescription, and other imple- 
ment a tion details (iCeverino et alJ l2010l; IGovernato et al.l 



20101: iKrumholz fe Dekell 120121: IScannapieco et al.1 l20l1 



Christensen et alJl2012l : iHopkins et al.ll2012h . 

Generically, we find that increasing the strength of feed- 
back results in the reduction of star formation efficiencies 
and the incre ase in galaxy sizes for a given stellar mass, as 
found in, e.g., ISales et al.1 (|2010h . Galactic outflows also in- 
crease the degree of rotational support of galaxies by prefer- 
entia lly removing low angular momentum gas from the ir cen- 
ters (IGovernato et al.ll2010l : iBrook et all 120111 . l2Q12h . Disk 
galaxies may survive or even be produced from high angular 
momentum mergers of gas-rich systems, provided that pres- 
surization from a multiphase ISM preven ts fragmentation 
and e fficien t conversion of gas int o stars ([Robertson et al.l 
120041 . 120061 : IHopkins et all l2009al lbh. Simulations with no 
thermally pressurized ISM, however, tend to produce more 
turbulent and clumpy disks and ma y result in rather differ - 
ent merger remnant morphologies ([Bournaud et al . 2011). 
These illustrate how the structure and morphology of high 
redshift galaxies provide yet another important test of feed- 
back processes during the peak epoch of galaxy growth. 

The progress of high performance computing has en- 
abled these types of cosmological zoom simulations to be- 
come more routine, but they still rely on subgrid prescrip- 
tions and are subject to numerical uncertainties. Advancing 
this field will rely on developing more robust models for the 
small-scale physics in addition to achieving higher dynamic 
range. This work presents a first step towards developing 
plausible models for the formation of galactic systems down 
to sub-kpc scales, but much work lies ahead to push both 
to smaller scales and to understand better the underlying 
physical processes governing high redshift galaxy assembly. 
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